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ABSTRACT 

We perform a statistical analysis of the temporal and spectral properties of the latest Fermi gamma-ray 
bursts (GRBs) to revisit the classification of GRBs. We find that the bimodalities of duration and the energy 
ratio (£peak/Fluence) and the anti-correlation between spectral hardness (hardness ratio (HR), peak energy and 
spectral index) and duration (Tgo) support the long/soft - short/hard classification scheme for Fermi GRBs. 
The HR - Tgo anti-correlation strongly depends upon the spectral shape of GRBs and energy bands, and the 
bursts with the curved spectra in the typical BATSE energy bands show a tighter anti-correlation than those with 
the power-law spectra in the typical BAT energy bands. This might explain why the HR - Tgo correlation is 
not evident for those GRB samples detected by instruments like Swift with a narrower/softer energy bandpass. 
We also analyze the intrinsic energy correlation for the GRBs with measured redshifts and well defined peak 
energies. The current sample suggests £p, res t = 2455 x (i?i SO /10 ) for short GRBs, significantly different 
from that for long GRBs. However, both the long and short GRBs comply with the same E^ K!it -Li so correlation. 

Subject headings: gamma-rays burst: general - methods: data analysis 



1. INTRODUCTION 

The field of gamma-ray bursts (GRBs) has rapidly ad- 
vanced in recent years especially after the launch of NASA 
missions Swift (in 2004, Gehrels et al. 2004) and Fermi (in 
2008, Atwood et al. 2009). A physical classification of GRBs 
is still a basic open question (e.g., Zhang 2011). According 
to the traditional classification schemes, GRBs can be divided 
into long and short ones, based on the well-known bimodal 
distribution of their durations monitored by the Burst And 
Transient Source Experiment (BATSE, Meegan et al. 1992), 
which also show different spectral hardness ratios (HR; Kou- 
veliotouet al. 1993). The HR in conjunction with the duration 
provides a means for classification, e.g., the long/soft class 
comprises roughly 3/4 of the population and the short/hard 
class comprises the other 1/4 for BATSE GRBs (Kouveliotou 
et al. 1993; Qin et al. 2000). The difference between long 
and short GRBs is further established by the observations of 
afterglows and host galaxies. The fact that several nearby 
long GRBs are associated with Type Ic Supernovae (SNe) 
and most long GRB host galaxies are found be dwarf star- 
forming galaxies favors the speculation that most long GRBs 
are accompanied by massive stellar explosions (see Woosley 
& Bloom 2006, for a review). Some nearby short GRBs (or 
short GRBs with a long-soft extended emission) have host 
galaxies that are elliptical or of early type, with little star for- 
mation (Fox et al. 2005; Gehrels et al. 2005; Berger et al. 
2005a; Barthelmy et al. 2005). This points towards a dif- 
ferent type of progenitor, e.g. compact object mergers (see 
Nakar 2007, for a review). 

This dichotomous picture was soon challenged by some 
following observations. GRB 060614 and GRB 060505 are 
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two nearby long GRBs that did not have bright SN associa- 
tions, sharing similar properties to short GRBs (Fynbo et al. 
2006; Delia Valle et al. 2006; Xu et al. 2009). Three high-z 
GRBs 080913, 090423 and 090429B have rest-frame dura- 
tions shorter than 1 s, but are likely related to massive star^l 
(Greiner et al. 2009; Salvaterra et al. 2009; Tanvir et al. 2009; 
Cucchiara et al. 201 1). An observed short GRB 090426 was 
found in many aspects similar to long GRBs and it was also 
probably linked to the death of a massive star (Antonelli et 
al. 2009; Levesque et al. 2010; Xin et al. 2011; Thone et al. 
2011; Nicuesa Guelbenzu et al. 2011). These results suggest 
that certain observational properties (e.g., long vs. short du- 
ration) do not always refer to certain types of progenitor (see 
Zhang et al. 2009 and the references therein). Moreover, the 
bimodal duration distribution was not presented in observa- 
tion for some detectors with a narrow/softer energy bandpass 
such as HETE-2 and Swift, and the distinction of short/hard 
and long/soft in the hardness-duration panel is not very clear 
(e.g., Sakamoto et al. 2011 and the references therein; Shao 
et al. 2011). Recently, Guiriec et al. (2010) analyzed 3 very 
bright short GRBs observed by the Fermi Gamma-Ray Burst 
Monitor (GBM, Meegan et al. 2009), and found that short 
GRBs are very similar to long ones, but with light curves con- 
tracted in time and with harder spectra stretched toward higher 
energies. They also showed that the hardness evolutions dur- 
ing the bursts follow their flux/intensity variations, similar to 
long bursts. By studying the composite light curves including 
both the prompt and afterglow emission of GRBs detected by 
Swift, Shao et al. (2011) found the similarity of the radia- 
tive features between the long and short bursts. They also 
proposed that the spectral evolution of the prompt emission 
might be an important factor that determines the correlation 
between the hardness ratio and duration. 

Theoretically speaking, the duration may be not an unique 
indicator for the physical category of a GRB. In the pre-Swift 
era, it was already known that (1) the coalescence of two com- 
pact objects can produce either long or short GRBs, depend- 

5 The other possibility that these high redshift bursts may be from the 
super-conducting cosmic strings has been tightly constrained (Wang et al. 
2011). 
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ing on the remnant formed in the merger (a short burst is likely 
if the remnant is a stellar black hole while a long burst is 
expected if the remnant is a super-massive neutron star; see 
Kluzniak & Ruderman 1998 and the references therein); (2) 
the collapse of a massive star usually produces long GRBs 
but the possibility of generating short events can not be ruled 
out (Zhang et al. 2003; Fan et al. 2005). One may appeal 
to multiple observational criteria to judge the correct physi- 
cal category of the GRB progenitor model that is associated 
with a certain GRB (e.g., Zhang et al. 2009; Fan et al. 2005; 
Kann et al. 2011). Recently, Goldstein et al. (2010) found 
that the distribution of the Epeak/Fluence energy ratio has a 
clear bimodality by analyzing the complete BATSE 5B Spec- 
tral Catalog. An obvious distinction between long and short 
bursts emerges from this bimodal distribution. This result was 
further confirmed by an analysis of 382 GRBs from the GBM 
spectral catalog (Goldstein et al. 201 1). Another phenomeno- 
logical classification method for GRBs was proposed by Lii 
et al. (2010). 

It is urgent to discuss the classification issues further in the 
Fermi era. The preliminary analysis results of GRBs observed 
by Fermi-GBM presented by e.g., Bissaldi et al. (201 1), Nava 
et al. (201 1, hereafter Nil), and Shao et al. (201 1) confirmed 
the duration bimodality found in the BATSE data (Kouve- 
liotou et al. 1993). Whether the hardness-duration correlation 
is also consistent with the BATSE result and the short/hard - 
long/soft dichotomy classification scheme is robust have not 
been explored yet. In this paper, we present a systematic anal- 
ysis of the temporal and spectral properties of Fermi-GBM 
GRBs catalogued by Nil and analyze the intrinsic energy 
correlation of GRBs with measured redshift and well defined 
peak energy. This paper is structured as follows. The data and 
sample are presented in Section 2. In Section 3, we revisit 
the distribution of duration, the correlation between hardness 
and duration, the distribution of energy ratio and the intrinsic 
energy correlation. In Section 4, we give a summary. 

2. DATA AND SAMPLE 

The Fermi-GBM detected 438 GRBs by the end of March 
2010. The spectral properties of these GRBs have been ana- 
lyzed and published by N 1 1 . Out of the 432 GRBs for which it 
was possible to perform the spectral analysis, 323 bursts have 
curved spectra which are well fitted by the Band (Band et al. 
1993) or by a cut-off power law (CPL) model. The remaining 
109 bursts are best fitted with a simple power-law. In addition, 
the peak flux spectra of 235 bursts could be extracted and fit- 
ted with a Band or CPL model. The detailed data extraction 
and analysis can be found in the catalog by Nil, which we 
adopted in the following analysis. Recently, the first two years 
Fermi-GBM GRB catalog was released (see, Paciesas et al. 
2012; Goldstein et al. 2012). The duration (r 90 ) of 487 GRBs 
are reported in the catalog. Combined with the duration data, 
we separate the 427 GRBs with an analyzed spectrum into dif- 
ferent samples (5 cases are not included for lack of data): 
322 GRBs with the curved spectra (Band or CPL model) rep- 
resent sample 7, while 103 GRBs with the power-law spectra 
represent sample Furthermore, the 234 GRBs with a fitted 
peak flux spectrum represent sample 3 (one burst without Tgo 
data is excluded). 

To analyze the intrinsic energy correlation of the different 
GRB classes, we collect the GRBs with known redshift and 

6 The two GRBs with unusually soft spectra (a = -3.52 for GRB 1001 12 
and a = -3.98 for GRB 100207) are excluded in the analysis. 



well defined peak energy up to the end of May 2011. This 
sample includes 110 long GRBs and 7 short GRBs detected by 
BeppoSAX, HETE-2, Swift, Suzaku and Fermi. Most of data 
are taken from Amati et al. (2008, 2009), Nava et al. (2008), 
Ghirlanda et al. (2009, 2010) and the references therein, as 
well as the GRB Coordinates Network (GCN). The isotropic 
gamma-ray energy (E lso ) and luminosity (L; so ) are calculated 
in the rest frame in the energy range 1 - 10000 keV. 

3. TEMPORAL AND SPECTRAL ANALYSIS 
3.1. Duration and Hardness 

Distribution of duration is crucial for the traditional GRB 
classification. We analyze the distribution of durations of 
Fermi GRBs in Nil catalog in detail. As shown in Figure 1, 
the duration distributions are bimodal either for all the bursts 
or for samples 1 and 2. The "short" and "long" GRBs are well 
separated in the log-normal plot, which is consistent with the 
previous findings (e.g., Bissaldi et al. 2011; Nil; Shao et 
al. 2011). For all the GRBs in Nil, we find a cental value 
p\ — -0.23 (i.e., Tgo ~ 0.59 s; with a standard deviation 
cr\ = 0.47) and p2 - 1-42 (i.e., Tgo ~ 26.3 s; with a standard 
deviation <T2 = 0.47) for short and long bursts, respectively. 
For sample 1 and sample 2, the best fits yield p.\ = -0.36 
(790 ~ 0.44 s; a-i = 0.46) and p 2 = 1.47 (T 90 ~ 29.5 s; 
cr 2 = 0.47) and p. x = -0.23 (T 90 ~ 0.59 s; a\ = 0.22) and 
P2 - 1-04 (Tgo ~ 11 s; cr 2 = 0.57), respectively. These 
results show that the duration bimodality is not affected by 
the spectral shape of GRBs, but the duration central value of 
long bursts in sample 2 is shifted to smaller value and the 
duration distribution covers a smaller range (see the middle 
panel of Figure 1). In other words, on average, the duration 
of long GRB with power-law spectrum is smaller than that 
of GRB with curved spectrum. It is well known that the du- 
ration strongly depends upon the sensitivity and the energy 
range of the instrument. The bimodal distribution of duration 
is either a robust result or a instrument selection effect. Us- 
ing the data of GRBs simultaneously detected by Swift-BAT 
and Fermi-GBM, Virgili et al. (201 1) found that the duration 
distributions in the BAT and GBM bands are all bimodal. It 
is essential to make a more detailed study explaining the fact 
that the bimodal duration distribution is not presented in ob- 
servations of HETE-2 and Swift. The future work is going 
to explore it. We also adopt the conventional Tgo — 2 s to 
separate short and long GRBs. 

The spectral hardness is an additional discriminator for the 
classification of GRBs. In this work, we firstly explore the 
correlation between HR and Tgo, where HR defines the flu- 
ence ratio between two broad energy bands. We simply di- 
vide the GBM energy band (10-1000 keV) into five energy 
bands, i.e. 10-25, 25-50, 50-100, 100-300, 300-1000 keV. 
Using the spectral data in N 1 1 , we are able to calculate HR in 
arbitrary two energy bands. We define four different HR mea- 
surements, namely HR1, between the 50-100 keV and the 25- 
50 keV energy bands, i.e. the typical BAT energy bands; HR2, 
between the 100-300 keV and the 50-100 keV energy bands, 
i.e. the typical BATSE energy bands; HR3, between the 300- 
1000 keV and the 10-300 keV energy bands; and HR4, be- 
tween the 100-1000 keV and the 25-100 keV energy bands. 

We find that the values of HR are significantly different in 
these energy bands. This is expected, since they strongly de- 
pend on the energy bands between which they are calculated. 
Figure 2 shows HR versus Tgo for 425 GRBs from Nil. From 
this figure, we find that there is an obvious tendency that short 
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GRBs have harder spectra than long GRBs. The correlations 
between HR and Tgo are reported in Table 1. We find that 
(1) for all 425 GRBs, HR and Tgo are all anti-correlated for 
all choices of HR (HR1-HR4), but the correlation coefficients 
as well as the slopes are different; (2) the values of HR are 
larger and the HR - Tgo anti-correlation is stronger in the typ- 
ical BATSE energy bands (the median value of HR is 2.21, 
the correlation coefficient is r — -0.41, the chance proba- 
bility is p = 2.8 x 1(T 18 , and the slope is b = -0.12) than 
those in the BAT energy bands (the median value of HR is 
1.56, r = -0.28, p = 5.2 x 10~ 9 , and b = -0.06); (3) if we 
consider short and long GRBs separately, the correlations are 
very weak or even negligible. 

In addition, the question arises whether the HR - Tgo anti- 
correlation of GRBs depends on spectral shapes. For a com- 
parison with the previous results from the BASTE and BAT 
observations, the correlations between HRl and T 9 q, and be- 
tween HR2 and Tgo for sample 1 and sample 2 are investi- 
gated. The results are shown in Figure 3 and Table 1. We 
find that the anti-correlation between HR and Tgo indeed de- 
pends on the spectral shape of GRBs, and the GRBs with the 
curved spectra have a more clear correlation than those with 
the power-law spectra in the same energy bands (for a de- 
tailed analysis, see Table 1). Now, the fact that there is no 
obvious correlation between HR and Tgo in the Swift-BAT 
sample as well as other samples detected by instruments with 
a narrow/softer energy bandpass such as HETE-2 (Sakamoto 
et al. 201 1) can be easily understood. This is because (1) the 
HR - Tgo anti-correlation depends upon the energy bands and 
the correlation in the typical BAT energy bands is weak; (2) 
the majority of spectra of BAT GRBs are best fitted by a sin- 
gle power-law model. This is due to the fact that BAT only 
covers a narrow energy band. The HR - Tgo anti-correlation 
for GRBs with the power-law spectra is not obvious; (3) for a 
given GRB sample, the values of HR are small and the values 
of Tgo are large for BAT-like softer detector, this makes the 
HR - Tgo anti-correlation unclear. 

The peak energy and spectral index are also usually adopted 
to depict the hardness of GRBs. Nl 1 found that, on average, 
short GRBs have a higher peak energy (~ 490 keV) and a 

harder low-energy spectral index ( 0.50) than long GRBs 

(~ 160 keV and ~ -0.92). We analyze the correlations be- 
tween £peak (peak energy in the time integrated spectra) and 
Tgo, and between (peak energy in the peak flux spectra) 
and Tgo for sample 1, and between a (which represents the 
spectral index of GRBs best fitted with a power-law model) 
and Tgo for sample 2. The results are shown in Figures 4 
and 5. We find that £ pea k and Tgo, £p eak and Tgo, a and Tgo 
are all anti-correlated. The detailed correlation analysis are 
presented in Table 2. The fact that short GRBs have harder 
spectra and long GRBs have softer spectra is confirmed once 
again. Recently, Gruber et al. (2011) suggested that the rest 
fame peak energy distributions might be the same for the two 
classes of GRBs based on a small Fermi GRB sample with 
measured redshift. The different redshift distribution of the 
two classes of GRBs might be responsible for the different 
distributions of the observed peak energies (see e.g, Guetta & 
Piran 2005). However, we also should note that the redshift 
distribution which has been found to be different for long and 
short GRBs might have been strongly affected by the mea- 
surement methods. Short GRBs tend to have lower redshift, 
very similar to those of long GRBs measured by the same 
method, i.e., spectral analysis of the presumed host galaxies 



(Shao et al. 2011). Therefore, the distribution of the intrin- 
sic peak energy of the two classes of GRBs and the £p ea k - 
Tgo anti-correlation are needed to a further confirm. More- 
over, for the short and long classes separately, the £ pea k - Tgo, 
E P eak - Tgo, and a - Tgo anti-correlations do not exist or even 
are positive correlated. These results are consistent with what 
was found from the analysis of the correlation between HR 
and Tgo. Our results indicate that hardness ratio, peak energy 
and spectral index might be linked to the same physical fea- 
ture of GRBs. The nature is unknown so far. 

3.2. Distribution of Energy Ratio 

Although the distributions of hardness and duration can be 
used to classify GRBs, the overlap between these two classes 
of GRBs can not be ignored. Moreover, this scheme strongly 
relies on the subjective choices required for the duration cal- 
culation. Recently, Goldstein et al. (2010) showed that the 
Zspeak/Fluence energy ratio (which physically represents a ra- 
tio of the energy at which most of the gamma rays are emit- 
ted to the total energy emitted in gamma rays) can be used 
as a new GRB classification discriminator. This also has the 
big advantage that it does not rely on the burst duration es- 
timate. Using the preliminary duration and spectral result of 
382 Fermi-GBM GRBs, Goldstein et al. (201 1) analyzed the 
distribution of energy ratio and found that the distribution sep- 
arated into long bursts and short bursts well. This supports 
the original claim obtained by Goldstein et al. (2010). Mean- 
while, we also analyze the distribution of the energy ratio for 
sample 1. The fluence in the 10-1000 keV energy band is 
calculated by using the spectral parameters in Nil, then we 
obtain the values of Ep^/Fluence for each of GRBs in sample 
1. Figure 6 shows the distribution of the energy ratio. In the 
top panel, we plot the 316 GRBs included in our sample 1. A 
bimodal distribution is evident, as previously shown by Gold- 
stein etal. (2011). By using a standard nonlinear least-squares 
fitting algorithm, we fit the distribution by two log-normal 
functions. Thebest fits yield a central value p\ = -1.34witha 
standard deviation o-\ = 0.62, and p2 = 0.27 with cr 2 = 0.35, 
respectively. In the bottom panel of Figure 6, we present two 
distributions corresponding to long (unfilled histogram) and 
short (filled histogram) GRBs. This further confirms the fact 
that (1) the bimodality of the energy ratio distribution is corre- 
lated to that of the burst duration, and that (2) the energy ratio 
is indeed a good discriminator for classifying GRBs. Indeed, 
it could be used to identify some of the controversial GBRs 
previously discussed (see Section 3.3). The energy ratio bi- 
modality can also be easily understood if we note that short 
GRBs tend to have larger peak energies and smaller fluences 
(due to their short durations) with respect to the long ones. 

To further check the differences between long and short 
GRBs, we investigate the correlations between the peak flux 
(P) and Tgo and between the E^JP ratio (which physically 

represents a ratio of the energy at which most of the gamma 
rays are emitted to the total energy emitted in gamma rays in 
one second at peak of a burst) and Tgo- Those two quantities 
namely do not depend on the duration of a burst. The detailed 
results are presented in Figure 7 and Table 2. Both the correla- 
tions between P and Tg and between £ peak /P and Tgo (shown 
in the top and bottom panel of Figure 7, respectively) are dif- 
ferent for long and short GRBs. An anti-correlation between 
P and Tgo and a positive correlation between £ peak /P and Tgo 
are found for all 234 GRBs with the peak flux curved spec- 
tra, although the correlations are not very strong. Likewise, if 
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we consider short and long classes separately, the correlations 
are negligible or even reversed. This further confirm that the 
correlation between the hardness and the duration of GRBs 
is only a general trend between two clusters of GRBs or two 
types of GRBs and does not apply to either type. 



3.3. Intrinsic Spectral Energy Correlation 

Amati et al. (2002) found a tight correlation between the 
rest frame peak energy, £p jres t and the isotropic equivalent 
gamma-ray energy, Ei so , which was confirmed by later obser- 
vations (Amati 2006, 2010). However, short GRBs do not fol- 
low the correlation, as is true for the peculiarly sub-energetic 
and close GRB 980425, the proto-type of the GRB/SN con- 
nection (e.g., Amati 2006, 2010; Piranomonte et al. 2008; 
Ghirlanda et al. 2009, Gruber et al. 2011). These facts sug- 
gest that the E V j^ t - E^o plane may be used to distinguish 
between different classes of GRBs and to understand the dif- 
ferences in the physics/geometry of their emission. Here we 
reanalyze this correlation taking into account the new obser- 
vational data. Only 7 short GRBs with redshift, reliable esti- 
mate of £peak and other spectral parameters are available by 
the end of May 201 1, as listed in Table 2. 

The three peculiar/controversial short GRBs (GRBs 
071227, 090927 and 100816A; Piranomonte et al. 2008; Am- 
ati et al. 2010; Gruber et al. 201 1) are excluded for the follow- 
ing reasons. Caito et al. (2010) suggested that GRB 071227 
represents another example of a disguised short GRB, after 
GRB 970228 and GRB 060614, on the basis of their analy- 
sis performed in the context of the fireshell scenario. GRB 
090927 and GRB 1008 16A are two short GRBs detected by 
Fermi and had been analyzed by Gruber et al. (201 1). How- 
ever the spectrum of GRB 090927 is adequately fitted by a 
simple power law function (see, Gruber et al. 2009; Nava et 
al. 2011), which has been confirmed by our current analy- 
sis. The category of this GRB hence can not be further iden- 
tified due to its inaccurate peak energy measurement. GRB 
1008 16A was simultaneously detected by Swift-BAT, Fermi- 
GBM and Konus-Wind, and its duration estimated by these 
three missions is 2.9+0.6 s (15-350 keV), 2 s (50-300 keV) 
and -2.8 s (20 keV-2 MeV), respectively. So it is difficult 
to determine whether this GRB belongs to the short or to 
the long class only based on its duration. Fan & Wei (2011) 
identified a possible wind-like medium surrounding this burst, 
which suggested a massive star origin (i.e, this short-like GRB 
should be a long one). Moreover, this GRB does not deviate 
from the E PtKSt - E^o region of long GRBs (see also Gru- 
ber et al. 2011). We calculate the energy ratio and obtain 
log(£ peak /Fluence)=-1.24 for GRB 100816A, which is a typi- 
cal value of long GRBs. Therefore, this event should belong 
to the long class and has been included in our long GRB sam- 
ple. 

We also collected GRBs observed by Fermi-GBM with the 
same standard as short GRBs analyzed above from the GCN. 
27 GRBs (including two short ones listed in Table 2) are ob- 
tained, as listed in Table 3. Using the same calculation method 
as Amati et al. (2008), we obtain their isotropic equivalent 
energies in the energy range 1-10000 keV (in the bursts' rest 
frame). Figure 8 shows the correlation between £p ;res t and 
Eis . Data reported in Amati et al. (2010) and the references 
therein are also included. We find that the short GRBs are 
significantly different from the long ones. To obtain a quanti- 
tative comparison, we fit the £p, rest - E liio correlations for the 



short and long GRBs separately. The best fits yield 

£ p , rest = 2455 x (-^)°- 59±004 for short GRBs, 

with a spearman's rank correlation coefficient r = 0.89 and a 
chance probability P = 6.8xl0~ 3 , and 

Ep , est = 100 x (^|)°- 51±0 03 for long GRBs, 

with r = 0.85 and P = 1.2xl0~ 31 (see Table 1). For long 
GRBs, the result is in agreement with that obtained by some 
authors (e.g., Amati 2010; Ghirlanda et al. 2009; 2010; Gru- 
ber etal. 2011). 

Wei & Gao (2003) discovered a tight correlation between 
the rest frame peak energy and the luminosity (L; so ) based on 
nine GRBs (see their Fig. 6). Such a correlation was soon con- 
firmed by others (e.g., Yonetoku et al. 2004; Ghirlanda et al. 
2005). Recently, further studies showed that although short 
GRBs are inconsistent with the £p a -est - Ei so correlation hold 
by long GRBs, they might follow the E PtKSt - Lj so correlation 
(e.g., Ghirlanda et al. 2009; Gruber et al. 201 1). Such a trend 
has been confirmed by our analysis with the latest datcQ (see 
Figure 9). The fit to the £p, re st - Li so correlation yields 

£ p , est = 302x(^)°- 4l)±003 , 

with r = 0.76 and P = 2.3xl0~ 23 . This might imply that the 
energy dissipation processes powering the prompt gamma-ray 
emission of short and long GRBs are rather similar though the 
progenitors of these two kinds of events are likely different. 

4. SUMMARY 

In this work, we perform a statistical analysis of the tem- 
poral and spectral properties of the latest Fermi-GBM GRBs 
catalogued by Nl 1 to revisit the classification of GRBs. The 
traditional short/hard and long/soft classification is supported 
for these Fermi GRBs by the following facts: (1) The duration 
of the prompt gamma-ray emission has a clear bimodal dis- 
tribution. (2) The energy ratio (i.e., £p ea k/Fluence) also has 
a clear bimodal distribution and its bimodality is correlated 
to that of the duration distribution. (3) The anti-correlation 
between spectral hardness (hardness ratio, peak energy, and 
spectral index) and duration is also confirmed. Moreover, we 
find out that the correlation between the hardness and the du- 
ration is only a general trend between two clusters of GRBs 
or two types of GRBs and does not apply to either type. In- 
terestingly there are two leading models for GRBs. One is the 
collapsar that likely produces long-lasting GRBs. The other is 
the merger of two compact objects which is expected to power 
short-living GRBs. However, it is not always reasonable to ar- 
gue a collapsar origin for long bursts and a merger origin for 
short bursts since the possibilities that the collapsar can also 
produce short GRBs and the merger can generate long GRBs 
can not be ruled out, as already realized in the pre-Swift era 
(e.g., Kluzniak & Ruderman 1998; Zhang et al. 2003; Fan 
et al. 2005). The observations of some Swift GRBs, such as 
GRB 060614 and GRB 080503 (e.g., Xu et al. 2009), partly 
confirm such speculations. Therefore additional discrimina- 
tors are highly needed to classify GRBs reliably. 

7 The long GRB data are taken form Ghirlanda et al. 2010 and the refer- 
ences therein. 
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In this work we also find that the anti-correlation between 
hardness ratio and duration depends upon the spectral shapes 
of GRBs and the energy bands. The bursts with the curved 
spectra in the typical BATSE energy bands (between the 100- 
300 keV and the 50-100 keV bands) show a tighter anti- 
correlation than those with the power-law spectra in the typi- 
cal BAT energy bands (between the 50-100 keV and the 25-50 
keV bands). This might explain why the HR - Tgg correlation 
is not evident for Swift-BAT GRB sample as well as other 
GRB sample detected by instruments with a narrow/softer en- 
ergy bandpass such as HETE-2 (e.g., Sakamoto et al. 2011, 
Shao et al. 2011). We also analyze the intrinsic energy cor- 
relation for the different GRB classes and find that all short 
GRBs deviate significantly from the £p res t - £)so correlation 



hold by long GRBs, and they might follow another one 



? p , res , = 2455 x (£ iso /10 3Z ) 



52^.0.59 



based on the current small sample. The future observation 
will test our result. 
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Fig. 1. — Distributions of the durations for 322 GRBs with the curved spectra (sample 1, top panel), 103 GRBs with the power-law spectra (sample 2, middle 
panel), and all 425 GRBs (bottom panel) observed by Fermi. The solid lines show the best fits with two log-normal functions and the dashed vertical line is 2 s 
separation line. The Tgo data are taken from Paciesas et al. 2012. 
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Fig. 2. — Hardness ratio versus duration for 425 Fermi GRBs. The filled circles represent the short GRBs and the filled triangles represent the long ones. The 
solid lines are the best fits for all bursts, the dashed lines are the best fits to short or long events separately. For clarity, the estimated errors are not shown. 
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Fig. 3. — The correlation between hardness ratio and duration for GRBs with the different spectral shape. The other symbols are the same as Figure 2. 
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Fig. 4. — Top: The correlation between the observed peak energies in the time integrated spectra (£ pea k) and durations for 322 Fermi GRBs with the curved 
spectra. Bottom: The correlation between the peak energies in the peak flux spectra (£p Cak ) and durations for 234 GRBs with the curved peak flux spectra. The 
other symbols are the same as Figure 2. 
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Fig. 5. — Distribution of the spectral indices (a) versus durations for 103 GRBs with the power-law spectra. The other symbols are the same as Figure 2. 
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Fig. 6. — Distribution of the energy ratio (£ pea k /Fluence), where the fiuence is measured in the 10-1000 keV energy range. Top panel shows the distribution for 
all 322 GRBs with the curved spectra. The dashed lines are the best fits with two log-normal functions separately and the solid line represent the superpositions 
of two log-normal fitting. Bottom panel shows the distribution for 275 long GRBs (unfilled histogram) and 47 short GRBs (filled histogram), where the solid 
curves are the best-fit log-normal functions. There are clearly two distinct distributions. 



12 



Zhang et al. 




Fig. 7. — The correlations between the peak flux (P) and the duration (Tgo), and between the E p ^/P ratio and Tyo for 234 GRBs with the curved peak flux 
spectra, where E 1 ' . is the peak energy in the peak flux spectra of a GRB. The other symbols are the same as Figure 2. 
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Fig. 8. — The correlations between the rest frame peak energy (£p, res t) and the isotropic total energy (Ei so ). The squares represent the short GRBs, the triangles 
represent Pre-Fermi long GRBs taken from Amati et al. (2010) and the references therein, and the circles represent Fermi long GRBs listed in Table 3. The stars 
are the two controversial GRBs, GRB 071227 and GRB 100816A. The solid lines are the best fit correlations: £ p-resl = 2455 X (-S^) " 59 for short GRBs and 

£ p ,re St = 100 X (-St) " 51 for lo ng GRBs. 
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Fig. 9. — The correlation between the rest frame peak energy (£p >res t) and the isotropic peak luminosity (Lj so ), where the 105 long GRB data are taken from 
Ghirlanda et al. 2010 and the references therein. The other symbols are the same as Figure 8. Combined short and long GRBs, the best fit to the £ pJ est - Li so 
correlation yield £ p , rcst = 302 X )0.40±0.03 
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TABLE 1 
Correlation analysis results. 



Data* 


Number 


HR* 


a 


b 


r 


Probability 


log(ffif) = a + b log(r 90 ) 


All 


425 


HRl 


0.27±0.01 


-0.06+0.01 


-0.28 


5.2X10-* 


All 


425 


HR2 


0.47±0.02 


-0.12+0.01 


-0.41 


2.8xl0~ 18 


All 


425 


HR3 


-0.15±0.05 


-0.29+0.04 


-0.42 


6.2X10- 20 


All 


425 


HR4 


0.68±0.03 


-0.25+0.02 


-0.44 


3.0xl0- 21 


Short 


77 


HRl 


0.30±0.02 


-0.08+0.04 


-0.22 


0.06 


Short 


77 


HRl 


0.52±0.02 


-0.09+0.05 


-0.13 


0.25 


Short 


77 


HR3 


0.03±0.06 


-0.09+0.13 


0.04 


0.74 


Short 


77 


HR4 


0.82±0.04 


-0.09+0.10 


-0.03 


0.81 


Long 


348 


HRl 


0.19±0.02 


-0.02+0.01 


-0.03 


0.56 


Long 


348 


HRl 


0.35±0.03 


-0.05+0.02 


-0.11 


0.03 


Long 


348 


HR3 


-0.43±0.11 


-0.11+0.07 


-0.14 


0.009 


Long 


348 


HR4 


0.46+0.06 


-0.10+0.04 


-0.15 


0.007 


SI All 


322 


HRl 


0.33+0.01 


-0.09+0.01 


-0.44 


8.4xl0" 17 


S2 All 


103 


HRl 


0.16±0.01 


-0.05+0.01 


-0.55 


1.5xl0- 9 


SI All 


322 


HRl 


0.49±0.02 


-0.14+0.02 


-0.36 


2.2x10-" 


S2 All 


103 


HRl 


0.42±0.01 


-0.07+0.01 


-0.55 


1.5xl0- 9 


SI Short 


47 


HRl 


0.39±0.02 


-0.02+0.04 


-0.09 


0.56 


S2 Short 


30 


HRl 


0.19±0.01 


0.02+0.03 


0.16 


0.40 


SI Short 


47 


HRl 


0.59±0.04 


-0.07+0.07 


0.05 


0.72 


S2 Short 


30 


HRl 


0.45±0.01 


0.02+0.03 


0.16 


0.40 


SI Long 


275 


HRl 


0.26±0.03 


-0.05+0.02 


-0.18 


0.003 


S2 Long 


73 


HRl 


0.13±0.02 


-0.03+0.02 


-0.20 


0.09 


SI Long 


275 


HRl 


0.32±0.04 


-0.03+0.03 


-0.07 


0.23 


S2 Long 


73 


HRl 


0.38+0.03 


-0.04+0.02 


-0.20 


0.09 



* SI represents sample 1 (316 GRBs with the curved spectra) and S2 represents sample 2 (108 GRBs with the power-law spectra). 

* The hardness ratios (HR) are measured between two different energy bands, namely HRl, between the 50-100 keV and the 25-50 keV energy 
bands, i.e. the typical BAT energy bands; HRl, between the 100-300 keV and the 50-100 keV energy bands, i.e. the typical BATSE energy 
bands; HR3, between the 300-1000 keV and the 10-300 keV energy bands; and HR4, between the 100-1000 keV and the 25-100 keV energy 
bands. 
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TABLE 2 
Correlation analysis results. 



Data* 


Number 


a 


b 


r 


Probability 


logCEpeak) = a + b log(r 9 o) 


SI All 


322 


2.59±0.03 


-0.21+0.03 


-0.30 


6.6xl0" 8 


SI Short 


47 


2.74+0.07 


0.01+0.13 


0.08 


0.59 


SI Long 


275 


2.34±0.08 


-0.05+0.05 


-0.05 


0.42 


log(£ p cak ) = a + b log(r 90 ) 


S3 All 


234 


2.64±0.04 


-0.17+0.03 


-0.22 


7.0xl0" 4 


S3 Short 


33 


2.83+0.08 


0.06+0.15 


0.14 


0.42 


S3 Long 


201 


2.37+0.09 


0.01+0.06 


0.02 


0.79 


log(-a) = a + b \og(T 90 ) 


S2 All 


103 


0.16+0.01 


0.05+0.01 


0.55 


1.5xl0" 9 


S2 Short 


30 


0.13+0.01 


-0.03+0.02 


-0.16 


0.40 


S2 Long 


73 


0.19+0.02 


0.02+0.02 


0.20 


0.09 


log(P) = a + b log(r 90 ) 


S3 All 


234 


-5.53+0.05 


-0.30+0.04 


-0.25 


9.0xl0" 5 


S3 Short 


33 


-5.29+0.12 


-0.15+0.23 


-0.09 


0.61 


S3 Long 


201 


-6.08+0.11 


0.05+0.07 


0.003 


0.97 


log(£ p cak /P) = a + b log(r 90 ) 


S3 All 


234 


-0.60+0.04 


0.14+0.03 


0.18 


0.006 


S3 Short 


33 


-0.66+0.08 


0.21+0.16 


0.10 


0.59 


S3 Long 


201 


-0.30+0.09 


-0.06+0.06 


0.01 


0.85 


log(£ p . rest )= a + b log(£ iso /10 52 ) 


Short 


7 


3.39+0.04 


0.59+0.04 


0.89 


6.8x10^ 


Long 


110 


2.00+0.01 


0.51+0.03 


0.85 


1.2xl0- 31 


log(£ p , res ,)= a + b logCZWICV 2 ) 




112 


2.48+0.03 


0.40+0.03 


0.76 


2.3xl0- 23 



* SI, S2 and S3 represent sample 1 (322 GRBs with the curved spectra), sample 2 (103 GRBs with the power-law spectra) and sample 3 (234 
GRBs with the curved peak flux spectra), respectively. 



TABLE 3 

Short GRBs with measured redshifts and spectral parameters. 



GRB 


z* 


a 


Peak flux 


Range 




-Ep.rest 


Fluence 


range 




Ref. 










(keV) 


10 52 erg/s 


(keV) 


(10- 6 ) 


(keV) 


10 52 erg 




050709 


0.16 


-0.53+0.12 


5.1+0.5E-6 


2-400 


0.05±0.01 


97.4+11.6 


0.4+0.04 


2-400 


0.0033+0.0001 


1 


051221 


0.5465 


-1.08+0.13 


4.6+1.3E-5 


20-2000 


6.42+0.56 


620+186 


3.2+0.9 


20-2000 


0.3+0.04 


1 


061006 


0.4377 


-0.62+0.2 


2.1e-5 


20-2000 


1.78+0.23 


955+267 


3.57 


20-2000 


0.2+0.03 


1 


070714 


0.92 


-0.86+0.1 


2.8+0.3 


100-1000 


1.4+0.1 


2150+1113 


3.7 


15-2000 


1.1+0.1 


1 


090510 


0.903 


-0.80+0.03 


80 


8-40000 


13.1+0.87 


8370+760 


30+2 


8-40000 


3.75+0.25 


2 


1001 17A 


0.92 


-0 14+ - 33 

-0 27 


6.1+0.4 


8-1000 


0.41+0.05 


551 + 142 


0.41+0.05 


8-1000 


0.09+0.01 


3 


101219A 


0.718 


-0 22 +0 - 30 


2.8+0.8E-6 


20-10000 


3.78±1.08 




3.6±0.5 


20-10000 


0.49+0.07 


4 


071227* 


0.383 


-0.7 


3.5±l.lE-6 


20-1300 


0.25±0.08 


1384+277 


1.6+0.2 


20-1300 


0.1 ±0.02 


5 



References: (1) Ghirlanda et al. (2008); (2) Guiriec et al. 2009 ; (3) Paciesas et al. 2010; (4) Golenetskii et al. 2010; (5) Golenetskii et al. 2007. 
Notes. * Some authors suggested that this GRB is a disguised short one as GRB 060614 (e.g., Caito et al. 2010). 
* The redshift are taken from the Greiner's webpage (http://www.mpe.mpg.de/~jcg/grbgen.html I. 
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TABLE 4 

The spectral parameters of 27 Fermi-GBM GRBs with known redshifts and well measured peak energies. The redshift are taken from the 

Greiner's webpage. 



GRB 



(keV) 



Fluence 
(1CT 6 erg/cm 2 ) 



range 



GCN 



(keV) 


(10 52 erg) 


number 


jU-jUU 


JJ. J±Z.+ 


o I uu 


25-1000 


2 27+0 76 


8263 


8-30000 


3 87 ±46 


8278 


25-900 


1 8+0 02 


8369 


8-1000 


28.8±1.7 


8715 


8-1000 


327±3 


9035 


8-1000 


19.3±0.2 


9057 


8-1000 


10.6±2.9 


9229 


8-1000 


4.47±0.09 


9230 


8-40000 


3.75±0.25 


9336 


s 1 nnn 


OS S + 1 Q 9 




8-1000 


25.4±0.6 


9535 


50-10000 


305±2 


9866 


8-1000 


186+5 


9933 


10-1000 


3.55±0.12 


9957 


8-1000 


10.6±0.1 


9983 


8-1000 


12.2±2.4 


10095 


8-1000 


1.63±0.02 


10204 


8-1000 


2.01 ±0.07 


10266 


8-1000 


0.09±0.01 


10345 


8-1000 


76.6±1.2 


10595 


8-1000 


2.66±0.11 


11015 


10-1000 


14.8±0.5 


11099 


10-1000 


0.73±0.02 


11124 


10-1000 


28.9±0.3 


11248 


10-1000 


0.59±0.04 


11477 


10-1000 


6.9±0.2 


11727 



080810 

0809 16A 

0809 16C 

081007 

081222 

090323 

090328 

090423 

090424 

090510" 

090516 

090618 

090902B 

090926A 

090926B 

09 1003 A 

091020 

091127 

091208B 

1001 17 A" 

100414A 

100728B 

100814A 

100816A* 

100906A 

101219B 

110213A 



3.35 

0.689 

4.35 

0.5295 

2.77 

3.57 

0.736 

8.26 

0.544 

0.903 

4.109 

0.54 

1.822 

2.1062 

1.24 

0.8969 

1.71 

0.49 

1.063 

0.92 

1.368 

2.106 

1.44 

0.8049 

1.727 

0.55 

1.46 



-0.91±0.12 

-0.9±0.1 
-0.91±0.02 
-1.4±0.4 
-0.55±0.07 
-0.89±0.03 
-0.93±0.02 
-0.77±0.35 
-0.9±0.02 
-0.8±0.03 
-1.51±0.11 
-1 26 +006 
-0.70±0.01 
-0.75±0.01 
-0.13±0.06 
-1.13+0.01 
-0.2±0.4 
-1.27±0.06 

-1 48+ - 05 
^ -0.05 

-0 14+ 033 
-0.58±0.01 
-0.9±0.1 

-0 64+ 014 
u -° -0.12 



-2.08±0.06 

-2.1±0.06 

-2.2±0.1 

-2.9±0.1 
-2.6±0.3 



9 5+0.15 
-0.33 
1 oc+0.21 

7 cq+0.04 
"^""-0.05 



-2.64±0.24 
-1.7±0.02 
-2.2±0.02 



9 O9+0.09 
~^- U -0.12 



-0.31±0.05 -2.77±0.17 



-1 34 +M 
• -0.06 

0.33±0.36 
-1.44±0.05 



-1 98+ 006 

i - :,o -0.07 
-2.12±0.12 



313.5±73.6 

109±9 

424±24 

40±10 

134±9 

697±51 

653±45 

82±15 

177±3 

4400±400 

185.6^ 

155 5 +nl 

-10.5 

775+11 

314±4 

91±2 

486.2±23.6 
47.9±7.1 
36±2 
144 2 +18 

287^ 
627.6! ||5 
131+15 ' 

106 4 +13 ' 9 
luu - -12.6 

136.7±4.7 

1UD -20.2 

70±8 

98 4 +8 5 
-6.9 



6.9±0.5 

15±5 

190 

1.2+0.1 

13.5+0.8 

100+1 

80.9+1 

1.1+0.3 

52+1 

30±2 

23±5 

270±6 

374±3 

145±4 

8.7±0.3 

37.6±0.4 

10±2 

18.7±0.2 

5.8±0.2 

0.41±0.05 

129±2 

2.4±0.1 

19.8±0.6 

3.84±0.13 

26.4±0.3 

5.5±0.4 

10.3±0.3 



Notes. " Two short GRBs. 

* This GRB was classified as a short one by the previous study (e.g., Gruber et al. 2011), we rind it should belong to the long class. 



